Copy 

RM L50F26a 




RESEARCH MEMORANDUM 


SOME CALCULATIONS OF THE LATERAL RESPONSE OF TWO 
AIRPLANES TO ATMOSPHERIC TURBULENCE WITH RELATION TO THE 

LATERAL SNAKING PROBLEM 
By John D. Bird 


Langley Aeronautical Laboratory 
Langley Air Force Base, Va. 



engineering dept, library 

CHANCE- VOUGHT AIRCRAFT 
DALLAS, TEXAS 


This document contains | 
affecting the national ! 
States within the meaning 
USC 50:31 and 32. 
revelation of its contenB 
unauthorized person isj 
Information so classi^ 
only to persons 
services of the Unitej 
civilian officers and e 
Government who hav« 
therein, and to United ^ 
loyalty and discretion * 
informed thereof. 


lified information 
of the United 
the Espionage Act, 
^ansmissi.in or the 
any manner to an 
•hibited by law. 
may be Imparted 
military and naval 
Jtates, appropriate 
pees of the Federal 
legitimate Interest 
Js citizens of known 
If necessity must be 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

WASHINGTON 

August 25, 1950 


NACA RM L50F26a 

NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


RESEARCH MEMORANDUM 

SOME CALCULATIONS OF THE LATERAL RESPONSE OF TWO 
AIRPLANES TO ATMOSPHERIC TURBULENCE WITH RELATION TO THE 
LATERAL SNAKING PROBLEM 
By John D. Bird 

SUMMARY 


Calculations are made of the lateral response to representative 
time histories of atmospheric turbulence for two airplanes having 
widely different dynamic properties, and explanations for their differ- 
ences in behavior are given. 

The results of the calculations indicate that under the proper con- 
ditions, atmospheric turbulence can initiate and maintain a lateral 
hunting oscillation of an airplane, and that this oscillation can be 
fairly regular in both amplitude and frequency. This effect is more 
pronounced for lightly damped airplanes. It is felt that this phenome- 
non may be the cause for some of the cases of airplane snaking that 
have not been explained by other considerations. 


INTRODUCTION 


Several recent high-speed airplanes exhibit a tendency to develop 
and maintain lateral hunting oscillations of roughly constant amplitude, 
which are generally referred to as snaking oscillations. These oscil- 
lations are of essentially the same period as the normal Dutch roll 
oscillation of the airplane and usually have amplitudes of the order 
of ±1° yaw. They are considered to be a nuisance to the operator, 
causing loss of confidence in -the airplane's response to control and 
making the airplane less satisfactory as a gun platform. 

Adequate explanations have been offered for this behavior in 
specific cases) however, there are still numerous occurrences for which 
no satisfactory explanation has been advanced. Some of the explanations 
for this motion are associated with nonlinear aerodynamic character- 
istics which result in different rates of damping for the large and 
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small amplitude ranges of motion. An example of this nature could 
arise from a poor fairing at the juncture of the tail surfaces. Other 
causes of snaking can be associated with small amounts of slack in the 
rudder control system, or the effects of fuel sloshing. 

Some unpublished experiments conducted in the Langley stability 
tunnel have suggested the possibility that motions, similar to those 
described as snaking, may result from the turbulence which exists in 
the atmosphere. The present paper, therefore, constitutes a preliminary 
investigation of this possibility. Calculations are made of the lateral 
response to representative time histories of atmospheric turbulence for . 
two airplanes, having widely different dynamic properties, and explana- 
tions for their differences in behavior are given. Calculations and 
experimental records of .the response of a model, which has freedom only 
in yaw, to the roughness in the stability-tunnel air stream are given 
for illustrative purposes. 


SYMBOLS 


The stability system of axes is used in the present analysis. This 
axis system has its origin at the center of gravity. The Z-axis is in 
the plane of symmetry and perpendicular to the relative wind, the X-axis 
is in the plane of symmetry and perpendicular to the Z-axis, and the 
Y-axis is perpendicular to the plane of symmetry. The positive direc- 
tions of the stability axes and of angular displacements of the airplane 
are ■ shown in figure 1. The coefficients of forces and moments employed 
in this paper are in standard NACA form. 

The coefficients and symbols are. defined as follows: 

Cp lift coefficient (Lift/qS) 

Cy lateral-force coefficient (Y/qS) 

C z rolling-moment coefficient (L/qSb) 

C n yawing -moment coefficient (N/qSb) 

dCv 

° Y P " 3T ' 

SC 7 

C *3 " W~ 
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C l n 


Sc n 

sr 


Sc 


z 


°2V 


Sc 


n 


-n T 


C Z, 


Cn T 


% 

K Z 

K XZ 

k* 

kz 

^b 

<P 

3 


SEk 
2 V 


SC 


Z 


A rb 

d 2V 


Sc 


n 

I 

2V 


s^ 


dimensionless radius of gyration about X-axis (k x /b) 
dimensionless radius of gyration about Z-axis (k z /b) 


product of inertia factor 


((* 


- k x I sin Ti cos T| 


radius of gyration about x-axis 
radius of gyration about z-axis 

W \ 


relative -density factor 


yy 


angle of bank of airplane, radians unless otherwise noted 
azimuth angle of airplane, radians unless otherwise noted 
angle of sideslip of airplane, radians unless otherwise 
noted (tan~ b - 


\ 


V 


angle of air stream with respect to initial flight -path 
. direction of airplane, radians unless otherwise noted 


amplitude of oscillation in air-stream direction 


k 
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*o 

1000 T|r 0 


amplitude of oscillation in model heading 
relative amplification, V in feet per second 


f 


t 

Dh 

s h 


frequency, cycles per unit time 
time 

differential operator ( d/ ds b ) 
dimensionless time (tV/h) 



d n 

ds b n 


w 

s 

b 

g 

X 

Y 

Z 

L 

N 

pb/2V 

rb/2V 

P 

r 


inclination of principal longitudinal axis of inertia with- 
respect to flight path} positive when the principal axis 
is above flight path at nose 

weight 

wing area 

wing span 

acceleration of gravity 

longitudinal force along X-axis 

lateral force along Y-axis 

normal force along Z-axis (Lift = - Z) 

rolling moment about X-axis 

yawing moment about Z-axis 

wing-tip helix angle, radians 

yawing-velocity parameter, radian measure 

rolling angular velocity about X-axis, radian measure 

yawing angular velocity about Z-axis, radian measure , 


v 


linear velocity of airplane along Y-axis 
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lateral gust velocity with respect to undisturbed position 
Of airplane 


free-stream velocity 
dynamic pressure 


mass density of air 

period of free lateral oscillation 

time for free lateral oscillation to damp to one-half 
amplitude 


CALCULATION METHODS 


General methods are given in references 1 and 2 for calculating 
the lateral response of aircraft to gusts. Reference 1 indicates that 
the response of an airplane to an arbitrary gust structure may be 
obtained by superposition of solutions for unit gusts, which, in the 
limiting case of a continuous disturbance function, involves the evalua- 
tion of Duhamel's integral. It is also pointed out in reference 1 that 
a rigorous analysis of gust effects requires consideration of penetration 
time and the aerodynamic lag in building up the lift on the surfaces. 
Examination of some of the penetration effects indicated that their 
magnitudes were small compared to the effects of sideslip. The rolling 
component of the measured turbulence considered for this paper was of 
course unknown, and thus all gust disturbances were necessarily con- 
sidered to be in a single plane. For the purpose of the calculations 
of this paper the turbulence was assumed to contribute nothing more than 
an effective change in sideslip of the airplane. Thus, in a side gust 

dC„ dC 7 

of velocity Vo, the angular and linear disturbances are o — — , a — 

do do 

dCy 

and o where the derivatives are numerically equal to C n g, C^g, 

° dCy 

and Cyp, respectively. In this analysis o — - was assumed to be zero. 


The lateral response of two airplanes to representative time 
histories of atmospheric turbulence was calculated by obtaining the 
motion of the airplanes following the application of unit yawing- and 
rolling -moment coefficients and then evaluating Duhamel's integral with 
this unit solution as the response variable and the record of the lateral 
fluctuation in air-stream direction as the forcing function'. 
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Now.Duhamel’s integral may be written 


^sbi = 


• s b^ 




(sbi - sb) 


da 

dsb 


dsb 


where ^a^b-^ - s-^ is the lateral response of the airplane to 
unit a. This equation may be broken into two integrals 




where 





dC Z da 
da ds^ 


ds> 


are the lateral responses of 


the airplane to unit yawing- and rolling-moment coefficients, respec- 
tively, and dC n /da and dC z /da are numerically equal to C n ^ 

and C 2 p , respectively. The last equation may be expressed as a single 
integral, of course 


^ s b x 





The solutions to the lateral equations of motion following the applica- 
tion of unit yawing- and rolling -moment coefficients were obtained by 
use of an automatic digital computing maching employing the procedures 
of reference 2. Duhamel^s integral was evaluated, by a numerical inte- 
gration, in a manner similar to reference 3 to obtain the motion of the 
airplane in response, to the turbulence. This calculation was also car- 
ried out on an automatic computing machine. 

A brief resume of the methods used for calculating the lateral 
frequency-response characteristics, of the airplanes considered is given 
in the appendix along with the equations of motion from which the 
response of the airplanes to unit disturbances was calculated. 


RESULTS AND DISCUSSION 
Calculated Airplane Motions in Turbulent Air 


The mass and aerodynamic characteristics of two airplanes having 
widely different operating conditions and dynamic characteristics are 
given in table I. Airplane A is a low-speed, low-altitude airplane and 
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airplane B is a high-speed, high-altitude research airplane, which is 
known to exhibit small continuous lateral oscillations under certain 
flight conditions. 

Figures 2 and 3 show the results of the calculation of the lateral 
response, of airplanes A and B to two known distributions of atmospheric 
turbulence. The turbulence shown in figure 3 was measured by. means of 
a sensitive recording pitot mounted on an airplane which traversed a 
region of turbulent air and so recorded the fluctuations in forward 
speed through the region. The turbulence shown in figure 2 was measured 
by means of an accelerometer mounted at the center of gravity of an air- 
plane which traversed a region of turbulent air. Both of these records 
of turbulence were considered to be fluctuations ' in sidewise velocity 
for these calculations. Reference 4 gives information on the measure- 
ment of atmospheric turbulence and justification for the assumption 
that the turbulence is isotropic. 

The gust distributions shown in the figures were assumed to exist 
in like fashion along the flight paths of both airplanes. The high- 
speed airplane, of course, encounters gusts with a greater frequency 
than the low-speed airplane. The gust velocities are assumed for these 
calculations to be the same for all altitudes. These results are thus 
of a qualitative nature. 

The distance traversed by the airplanes is used as an abscissa in 
the plots given, and the azimuth angles of the airplanes are chosen, to 
indicate the oscillation performed. The azimuth angle should be roughly 
proportional to the, apparent lateral' movement of the horizon. 

The motions of the two airplanes in response to the two regions of 
atmospheric turbulence are markedly- different (figs. 2 and 3). Air- 
plane A shows a response which might logically be termed by the operator 
as rough air; that is, the airplane responds in almost direct proportion 
to the local gustiness and subsides to little or no motion as the gusti- 
ness subsides. The response of most aircraft in the past seems to have 
been of this nature. Airplane B shows a response which builds up to a 
fairly steady lateral oscillation, which is almost independent of the 
local turbulence. This motion is very similar in character to motions 
which have been termed "snaking." 

The periods of the lateral oscillations, shown in figures 2 and 3, 
are close to the period of the classical free -lateral oscillations of 
the airplanes (table I). This is more nearly the case for airplane B 
than airplane A. 

Figure 4 shows the motions of both airplanes compared with the 
angular motion of the air with respect to the undisturbed attitude of 
the airplane for the case of the turbulence determined from airspeed 
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fluctuations. It should be noted that the turbulence record corresponds 
to different amplitudes of a for the two airplanes because of differ- 
ences in forward speed. It is easy to note a closer approximation in 
the case of airplane A to a one-to-one correspondence of air motion to 
airplane motion than for the case of airplane B, where there is little 
apparent relation between the air and airplane motions. It should be 
mentioned here that the long period change in heading shown for air- 
plane A in figure 3 was modified to some degree in the preparation of 
figure 4 in order to have the air and airplane motions oscillate about 
the same mean and so make for an easier comparison of the two motions. 

A comparison of a part of the lateral motion of airplane B, as 
calculated from the atmospheric turbulence (fig. 3)> with the snaking 
of this airplane recorded during a flight test is shown in figure 5. 

This comparison merely confirms the statement made previously that 
lateral motions arising from turbulence in the air can be similar in 
nature to flight measurements of snaking motions. The flight conditions 
for the two motions given are not identical, and the atmospheric turbu- 
lence existing during the flight test is unknown. The indications are, 
however, that turbulence having about one -third the magnitude of that 
employed for figure 3 would be required to maintain a hunting motion of 
airplane B comparable to the snaking. motion that it exhibited in flight. 
It should be pointed out that the disturbances do not have to be of the 
type generally referred to as "sharp-edged" gusts as shown in figure 2, 
but may be of a more gentle nature, as shown in figure 3- 


Free Motions of Model in Wind Tunnel 

As an example of the response of a free body to turbulence in the 
air stream, figure 6 gives the experimental and calculated hunting 
motion of a model mounted with freedom only in yaw in the air stream of 
the Langley stability-tunnel test section. The measured time history 
of the air-stream azimuth angle from which the model motion wqs calcu- 
lated is also shown. The calculation procedure was much the same as has 
been given previously for airplanes A and B. The time history of air- 
stream direction was obtained by use of a recording electronic pitot and, 
although the percentage error in the magnitude of the air- stream angles 
may be fairly large, the nature of the fluctuations should be accurate. 
The experimental model motion was not recorded at the same time as the 
air-stream azimuth angle. This fact does not change the nature of the 
result, however, because the air-stream turbulence was of a similar 
nature for a number of different recordings. 

The tunnel model was free only in yaw, was mounted on flexure plates 
in order to minimize friction, and consisted of only a fuselage and 
vertical tail. The mass and dimensional characteristics of the model 
oscillating system are given in table II. 
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It can be seen that the experimental and calculated motions of 
figure 6 are of a s imi lar nature, indicating that the air-stream turbu- 
lence is a significant factor in the hunting motion experienced by the 
model. • 


Frequency-Re sponse Character i st ic s 

The response of airplanes A and B, and of the tunnel model mounted 
with freedom in yaw, to sinusoidal forcing functions of various fre- 
quencies is given in figure 7 • The forcing functions are in the form 
of changes in heading of the approaching air stream. The results are 
given in terms of the amplitude of motion in radians induced by a sinus- 
oidal lateral gust distribution having an amplitude of 1000 feet per 
second and in terms of the phase lag of the motion behind the forcing 
function. 

An examination of these curves indicates the source of the differ- 
ence in response of airplanes A and B to atmospheric turbulence. Air- 
plane B has a frequency response- not unlike the characteristics of an 
electronic band-pass filter which excludes those harmonics of the 
applied frequency which are very much different from the resonant fre- 
quency. A good measure of the selectivity of response of these airplanes 
is the ratio of the amplification at the natural frequency to the ampli- 
fication at very low frequencies. Airplane A responds to a greater 
degree than airplane B to those frequencies that are different from the 
resonant frequency; thus, the tendency for the one-to-one correspondence 
of air direction to airplane azimuth angle shown for airplane A in fig- 
ure k* In general, the sharper .the frequency-response curve the more 
nearly the response to atmospheric turbulence approaches a sinusoidal 
motion. ^ 

Any effect which reduces the rate of free damping of an airplane 
should tend to increase the peak of the frequency-response curve and 
make the phase -angle shift at the natural frequency more abrupt. Air- 
plane B, of course, has a low rate of damping for the condition inves- 
tigated herein (table I). 

The rate of free damping can be affected to a marked degree by a 
change in one or more of the aerodynamic stability derivatives. This 
fact would make calculations, based on estimated derivatives, of a 
questionable nature unless some form of check is available. The . sta- 
bility derivatives used for the calculations of this paper are believed 
to be reasonably accurate because the experimental and calculated rates 
of and P of the, lateral oscillation compare well (table I). 

The frequency-response curve of the oscillating model mounted in 
the air stream of the Langley stability tunnel is between the curves of 
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airplanes A and B with regard to selectivity but large compared with 
both as regards over-all response (fig. 7) • 


CONCLUDING REMARKS 


• The results of the calculations indicated that under the proper 
conditions, atmospheric turbulence can initiate and maintain a lateral 
hunting oscillation of an airplane, and that this oscillation can be 
fairly regular in both amplitude and frequency. This effect is more 
pronounced for lightly damped airplanes. It is suggested that this 
phenomenon may be the cause for some of the cases of airplane snaking 
that have not been explained by other considerations. 

Langley Aeronautical Laboratory 

National Advisory ..Committee for Aeronautics 
Langley Air Force Base, Va. 
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APPENDIX 


CALCULATION OF LATERAL FREQUENCY RESPONSE 

The lateral frequency response of an airplane to an imposed sinus- 
oidal variation of wind direction may be calculated by solving the 
standard lateral equations of motion of an airplane (with the proper 
forcing terms added) for the- steady- state motion. For a unit sinusoidal 
variation in air-stream direction the amplitude of this motion becomes 
the amplification factor or the amount that the air-stream fluctuation 
is magnified. ' o 

Within the limits of the approximations discussed, the lateral 
equations of motion of an airplane experiencing a sinusoidal variation 
in air-stream direction are 

/ 2 2 1 \ / PI \ fb 

(^M-b^x p ^Zp-^bJ^ + ^Pb^xz^b “ 2 ^Zr^b/ “ *"Zp^ = °o^Z (3 s ^ n ~v~ s b 

( 2 ^b K XZ D b 2 " | C n p D b) ^ + ( 2 » i b K Z 2l> b 2 " \ c n r D b) 1/ “ C n p P = a o C rip sin 2lt ^ s b 

- C L<P + 2 ^bV + - C Y p ) 3 =' 0 

where the terms C^p and C^p on the right side of the first two equa- 
tions are considered to be the equivalents of and — - 

'da da ' 

. The variation of side force with rolling and yawing velocity, a 
term associated with the glide-path angle, and the side-force forcing 
function are omitted in these equations. Calculations showed these 
factors to be of little importance. Replacing the two right-hand terms 
of these equations by unity gives the equations from which the unit 
solutions were obtained for use in the calculations of response to 
arbitrary turbulence . 

Solving the equations for the steady- state motion gives for the' 
azimuth angle \J» 
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where e is the angle of lag of the motion behind the disturbance, and 

for o 0 = 1.0, the term \fa ?■ + b 2 is the. amplification factor for the 
given imposed frequency or 

*2 » \Ja 2 + b 2 
a 0 

Now the amplitude of the yawing motion for a given amplitude of lateral 
gust velocity is 



This term is called the relative amplification when 1000 is substituted 
for Vgyg^. and V is given in feet per second. Now 


ke^e “ k-u^u 



and 


b = 


keC u + k u Ce 

p 2 + p 2 


where 

k e = nigfg 2 + m 0 

= m 3 f s 3 + m l f s 
C e = Bf s 4 - Df s 2 

e u “ ' M s 5 + cf s 3 - Ef s 


• and 


' _ 2jtfb 
s v 




0 
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m2 = 


m l = (I %%% - \ C np C Yp C l 3 ) a o 

K + 2u b k x 2c Y^C np - ^b c n p + 2 ^b k XZ G Yp) C z'p| cr o 
m 3 ^b 2k X 2c np + ^b 2k XZ C Zp) a o 


and where the expressions for A } B, C, D, and E are given in 
reference 5* 


The lag angle is given by 

.-1 a 


6 = COS 


i 


a 2 + b 2 


= sin 


-1 k 


l 


+ b^ 


Substitution of various values of the imposed frequency in the above 
expressions gives the frequency-response curve. 
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TABLE I 

MASS AMD AERODYNAMIC CHARACTERISTICS OF AIRPLANES 





Airplane A 

Airplane 

Pb 

• 

• . • 

16.8 

106.3 

% 2 



0.0061 

0.0051 

h 

• • • 


0.0264 

0.0409 

K xz • • 

• « • 


0.0058 

-0.0006 

t>, feet 

• • « 


33.6 

28 

Weight, pounds 



8,700 , 

11,050 

S, square feet 

• • # 

« • • 

250 

130 

C L * • • • • 

• • • 

• • • 

0.551 

0.343 

c 7 

• « • 

• • • 

-0.280 

-O.474 

c n 

n p 


• • • 

-0.085 

0 

C 2 r 

• • • 

• • ft 

0.090 

0.224 

S • • • 

... 

ft ft ft 

. . -0.270 

-0.170 

C lp ’ 


• mm 

- 0.049 

-0.101 

C np 


• ft • 

0.097 

0.217 

c Yp 

• • • 

• • • 

.. -0.665 

-O.878 

V, feet per second 

• • « 

• • • 

257 

746 

Altitude, feet ... 

• . • 

• ♦ • 

7,500 . 

30,000 

P (calculated), seconds . . . 

• • • 

• • • 

3.7 

1.48 

P (flight), seconds 

• • • 

• ft • 

3.5 

1.80 

T l/2 (calculated), seconds . . 


• • « 

2.5 

6.56 

T l/2 (flight), seconds .... 

• • . 

• • • 

3.0 

6.50 














y/ 

1 ^ 
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TABLE II 

MASS AND AERODYNAMIC CHARACTERISTICS OF MODEL 
MOUNTED WITH FREEDOM IN YAW IN 
STABILITY-TUNNEL AIR STREAM 


p 

IzZ-? foo ‘ t pound second 
q, pound per square foot 

hi , feet ....... 

V, feet per second . . 

■*-S, square feet .... 
C n ^, per radian . . . . 


0.3 

bO 

3 

186 

2.25 

0.063 


Cn r ;..... -0.13 

P (calculated), seconds . . . . 0.8 

T 1/2 (calculated), seconds b.b 


and S are given as dimensions upon which aerodynamic coefficients 
are based for this model. 
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Relative wind, V 






disturbance 
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Figure A.- Comparison of lateral motions of airplanes A and B with 
fluctuation. in air-stream direction. Turbulence distribution of 
figure 3. 


Azimuth angle Angle 0 f sideslip 

ft, degrees /3 , degrees 
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(b) Calculated response of airplane to atmospheric 
turbulence (from figure 3) 

Mach no. 0.74, altitude 30,000 feet, weight 1 1, 000 pounds 

Figure 5.- Comparison of calculated hunting motion of airplane B in 
response to atmospheric turbulence with measured hunting motion of 
airplane B obtained from a flight test. 
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Figure 7.- Lateral- frequency response characteristics of airplanes 
and B and model mounted with one degree of freedom in stahility- 
tunnel air stream. 








Airplane 



(b) Phase lag of motion behind impressed moment. 




